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ABSTRACT 


This thesis considers the problem of evaluating a producer's program in statistical 


process control, from the standpoint of the consumer. A model is postulated reflecting 


the variability in proportion nonconforming of a process and the characteristics of the 


final control chart in the process. From this the steady-state solution to a Markov chain 


is used to find the output proportion nonconforming of the process. 
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THESIS DISCLAIMER 


The reader is cautioned that computer programs developed in this research may not 
have been exercised for all cases of interest. While every effort has been made, within 
the time available, to ensure that the programs are free of computational and logic er- 
rors, they cannot be considered validated. Any application of these programs without 
additional verification is at the risk of the user. 
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I. «sNTRODUCTION 


The commander of the Naval Weapons Support Center in Crane, Indiana is con- 
cerned with the quality level of projectiles stored at this base. A projected manpower 
shortage there may require the use of alternative methods to assure the quality of these 
projectiles with fewer personnel. The base commander also questions the capability of 
the current methods being employed to provide him with accurate estimates of quality. 

For many years, acceptance sampling has been employed as a primary tool to pro- 
vide assurance that high quality of goods is accepted by the military. Some military 
standards provide the sampling procedures, sample size, and decision rules by which a 
given shipment of goods, usually called a lot, is either accepted or rejected. These 
standards typically require that the sample or samples be drawn randomly from the lot 
of goods. These methods have several weaknesses. One drawback is that an entire lot 
must have been produced before a decision can be made to accept or reject it. Ifa lot 
is rejected, delays may occur before acceptable goods can be acquired. In some cases 
the producer’s own quality programs are so effective that acceptance sampling may be 
unnecessary due to such a high rate of acceptance. One alternative to acceptance sam- 
pling is to conduct a program that reviews and evaluates the quality control procedures 
already in place at the manufacturer's facility. If the quality personnel at the Navy base 
approve the procedures used by the manufacturer, then they may choose to accept the 
goods from the manufacturer’s opcration without further inspection. In this thesis we 
examine one method for evaluating a quality control procedure. Approval of a pro- 
ducer’s quality control procedures could be as effective iis acceptance sampling in pro- 
viding quality to the Navy and might reduce the number of quality assurance personnel 
needed at the Naval Weapons Support Center. 

In Chapter I] we examine the current quality methods that are being used at the 
base for acceptance sampling. 

In Chapter III, we discuss the differences between statistical quality control and 
statistical process control. This leads us to the exploration of using control charts as a 
means of helping to control the level of quality in a process. To examine the effect of 
a control chart on a typical manufacturing process, we first create a stochastic model 


of a manufacturing process. The model employs a Markov chain to describe how the 


behavior of the process is affected when a control chart procedure is imposed on it. The 











decision rules of the control chart are embedded in the model. A small computer pro- 
gram encodes the model. Given the necessarv inputs, one can use this program to 
compute the output quality level of the manufacturing process. The results of the pro- 
gram are included in this thesis. They provide a rough measure of how the process and 
control chart parameters might impact the quality of the output of the manufaciuring 
process. 

In Chapter IV. cost and effectiveness measures are used to determine whether or not 
the proposed method of statistical process control could be more desirable than the 
<urrent method of statistical quality control. The cost comparison is essentially a qual- 
itative argument. The effectiveness comparison is quantitative, 

In Chapter V the results of the study are summarized and recommendations are 
made pertaining to the use of alternative quality assurance methods at the Naval 
Weapons Support Center. 
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Il, ACCEPTANCE SAMPLING: DESCRIPTION AND ANALYSIS 


In this chapter, the importance of an effective quality assurance program is ad- 
dressed. Also, current methods used in acceptance sampling are explored with the pri- 
mary focus being on MIL-STD 105D. MIL-STD 414, another set of acceptance 
sampling plans, is also discussed and comparisons are made between MIL-STD 105D 
and MIL-STD 414. 

Acceptance sampling may be used by a consumer in deciding to accept or reject a 
shipment of goods from a manufacturer. If the shipment of goods is rejected, it will of- 
ten be sent back to the manufacturer. If a manufacturer's product is being rejected at 
a high rate, one of two things may happen. The manufacturer may take steps to im- 
prove his production methods, or the customer may find a better source of supply. Be- 
sides helping to ensure quality in accepted goods, acceptance sampling indirectly 
improves quality of production through its encouragement of good quality by requiring 
a high rate of acceptance. 


A. MIL-STD 105D SAMPLING PLANS 

Some of the most popular acceptance sampling plans are those in MIL-STD 105D. 
MIL-STD 103D was developed in 1950. Many of the sampling tables and procedures 
used in MIL-STD 105 were derived from the Army Service Forces Tables developed in 
1942 by Bell Telephone Laboratories for the Ordnance Department of the United States 
Army. Another source for the development of MIL-STD 105 was the JAN (Joint Army 
and Navy) Standard 105 Statistical Sampling Tables and Procedures originally developed 
for the Navy by the Statistical Research Group of Columbia University in 1945. 
MIL-STD 105D. the fourth revision to the standard, was adopted in 1963 by the ABC 
Working Group, a committee made up of members from the militar agencies of the 
United States, Great Britain. and Canada. This was the last revision to the standard. 
[Ref. 1] 

MIL-STD 105D 1s essentially a collection of sampling plans or a sampling scheme. 
The plans in this standard only use attributes data. The focal point of the standard is 
the acceptable quality level or AQL. The AQL is “the maximum percent defective (or 
the maximum number of defects per hundred units) that. for purposes of sampling in- 
spection. caa be considered satisfactory as a process average “ (Ref. 2]. In applying the 


standard. it is expected that there will exist a clear understanding between the producer 











and the customer as to what the customer considers the acceptable quality level for a 
given product. Although an obvious purpose of an acceptance sampling plan is to’en- 
sure that the customer receives a product of at least acceptable quality, an effect of using 
the plan “is in general to force the supplier to submit product of such a quality that a 
small percentage of the lots submitted for inspection are rejected” [Ref. 1]. 

Another characteristic of the standard is inspection level. The inspection level de- 
termines the relationship between lot size and sample size. Three general levels of in- 
spection are offered. Under ordinary circumstances, Level II is used. In instances, 
where the quality of goods is high, we may specify Level I, which will provide less dis- 
crimination via a smaller sample size. When the quality of goods is low, we may specify 
Level III, which will provide more discrimination via a larger sample size. There are also 
four special inspection levels offered: S1, $2, $3, and $4, which offer smaller sample sizes 
than Level I. These special levels may be used where relatively small sample sizes are 
necessary and large sampling risks can or must be tolerated. The inspection level is 
adopted at the initiation of the sampling program and is generally not changed there- 
after. [Ref. 2] 

MIL-STD 105D also offers three different types of sampling plans to choose from: 
single sampling, double sampling, and multiple sampling plans. The choice of a plan is 
frequently based on cost comparisons between the difficulties in conducting the plans 
and the average sample sizes of the plans. ‘The average sample size for single plans is 
generally greater than that of either double or muluple plans. However, the adminis- 
trative difficulty and average cost per unit of the sample are usually less with the single 
plan than with either the double or multiple plans. [Ref 1! We only consider single 
sampling plans in this thesis. 

Vora specified AQL. inspection level, and lot size. MIL- STD 105D gives a normal 
sampling plan that is to be used as long as the supplier is producing a product of AQL 
quality or better. It also gives a tightened plan to shift to if there is evidence that a de- 
terioration in quality has occurred. The rule is that a switch from the normal plan to 
the tightened plan will be made if two out of five consecutive lots have been rejected on 
original inspection. Normal inspection is re-instituted if five consecutive lots have been 
accepted on original insrection. If ten consecutive lots remain under a tightened plan, 
inspection is stopped pending action on quality. [Ref 1] Similarily, MIL-STD 105D 
culls for shifts to reduced inspection if the quality is observed to be especially good. The 
rule here is that production must have been running at a steady rate and the last 10 lots 


must have been accepted on original normal inspection. In order to shift to reduced 
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inspection, the cumulative number of defective items in the last 10 samples must be less 
than the value sct forth in the MIL-STD 105D table pertaining to limit numbers for re- 
duced inspection, [Ref. 1] The basic idea of reduced inspection is to save the consumer 
money when quality has been consistently high. 

The following is an example of one sampling plan from MIL-STD 105D. Suppose 
lots of size 1000 have been specified, an AQL of 1.5 percent has been agreed to, and a 
discrimination Level II is desired. For lot sizes of 1000, the single sampling plans set 
forth in MIL-STD 105D require a sample of size 80. The lot is accepted if there are 3 
or fewer defective items in the sample and rejected if there are 4 or more defective items 
[Ref. 3]. Each sampling plan in MIL-STD 105D generates an operating characteristic 
(OC) curve. An OC curve gives the probability of accepting a lot of a given quality. 
Just as with any OC curve, there are Type ] and Type I] errors involved. In our case, 
the probability of a Type I error is the probability of rejecting a lot as a function of the 
quality of the lot (proportion non-conforming). The probability of a Tvpe II error is the 
probability of accepting a bad lot. where “bad” corresponds to a larger proportion 
nonconforming that the specified AQL. The OC curve for our example is shown in 
Figure 1]. 
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Figure 1. Acceptance Sampling OC Curve 


For this OC curve, « , or the probability of rejecting a lot of good quality (AQL 
= 015) is equal to approximately .01. Similarly, B , or the probability of accepting a 
bad lot(p = .06, for example) is equal to approximately .25, 


B. MIL-STD 414 SAMPLING PLANS 

Another frequently used set of sampling plans is found in MIL-STD 414. All of the 
sampling plans in this standard use variables data, and it is assumed that the variables 
have a normal distribution. It is one of several variables sampling plans. This standard 
is not used at the Naval Weapons Support Center. MIL-STD 414 is an example of a 


variables sampling plan for situations where the quality requirement is still specified as 
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specified as a percent nonconforming. [Ref 4] 
Some points of similarity between MIL-STD 105D and MIL-STD 414 are: 
® both are based on the concept of AQL, 
¢ both pertain to lot-by-lot acceptance inspection, 
® both provide procedures for normal, tightened, and reduced inspection, 
® both provide several inspection levels, 
¢ in both sets of plans, the sample size is greatly influenced by lot size, and 


¢ in both sets of plans, all plans are identified by a sample size code letter. 


MIL-STD 414 is also similar to MIL-STD 105D in that in both sets of plans sampling 
plans under normal inspection were designed to protect the producer by rejecting with 
small probability those lots which have a percent defective no larger than the specified 
AQL. Just as in MIL-STD 105D, the protection to the customer depends largely on the 
use of tightened inspection whenever the process average appears to be unsatisfactory 
or there are other reasons to be suspicious of the process. [Ref. 3] 

The primary advantage of using MIL-STD 414 ts that the same OC curve can be 
obtained with a smaller sample size than is required by Standard 105D. The precise 
measurements required by a plan from Standard 414 will probably cost more than the 
simple classification of items required by a plan from MIL-STD 105D, but the reduction 
in sample size will hopefully offset this extra expense. For example, a plan from 
MIL-STD 105D may require a sample of 50, but a comparable plan from Standard 414 
requires a sample of only 30. If the unit cost of measurement is less than 5 3 times that 
of simply classifving the items, the MIL-STD 414 sampling plan will effect a savings. 
[Ref. 3] Such savings may be especially significant if inspection is destructive and the 
item is expensive. 

The primary disadvantage of a sampling plan from MIL-STD 414 is that a separate 
plan must be employed for each quality characteristic that is being inspected. For ex- 
ample, if an item is to be inspected for five quality characteristics, then it would be nec- 
essary to have five separate inspection plans from MIL-STD 414: whereas, acceptance 
or rejection of the lot as a whole could be based on a single sampling plan from Standard 
105D. Also, it is theoretically possible, although unlikely, that under a variables plan a 
lot will be rejected by the variables criteria, even though the sample actually contains 
no defective items. Another disadvantage of a variables sampling plan is that the quality 


characteristic must come from a normal distribution. Since the quality characteristics 











being sampled will not always be normally distributed, the use of MIL-STD 414 will not 
always be appropriate. 

Due to the higher degree of popularity and use of MIL-STD 105D as opposed to 
MIL-STD 414, plus the fact that Standard 105D is used at the Naval Weapons Support 
Center, MIL-STD 105D will be used as the benchmark with which to compare any 
proposed methods in later chapters. 


C. MEASURING THE PERFORMANCE OF A SAMPLING PLAN 

The most basic measure of a sampling plan’s performance in monitoring quality is 
its OC curve, which shows the probability that a lot of given quality will be accepted. 
It is also useful to have a measure of the quality of the product that will reach the con- 
sumer’s shelves after acceptance sampling has taken place. 

We will assume that rectifying inspection is employe 1 the sampling process. 
What is meant by this is that all rejected lots will be inspect sercent and noncon- 
forming items will be replaced by good ones. Thus, with th’ imption, as lots of size 

L_ pass through the sampling process, one of the following scenarios will occur: 

1. The lot is accepted, or 

2. The lot is initially rejected; but. is subjected to 100 percent inspection and bad items 
are replaced with good ones. The lot is then accepted. 

The sample s will eventually contain no defective items since all defects in the 
sample are replaced by good ones. The same is true for the lot if it is initially rejected. 
It will then be inspected 100 percent with the bad items being replaced by good ones. 
If the lot is initially accepted, it will contain p(Z—s) defective items, where p is the 
proportion of defective items in the lot. This will occur Pa proportion of the time. 
Thus, the sampling plan will on the average, turn out lots that contain Patp(L —s)) 
defective items (Pa is the probability of accepting a lot and is taken from the OC curve 
for the sampling plan). If we divide by the lot size, Z we can express this average out- 
going quality as a fraction defective as follows: 


AOQ = [Ref 1]. 


Pa(p(L — s)) 
L 
Thus, this concept of average outgoing quality, or AOQ, gives us a method of examining 
the effect that a particular acceptance sampling plan has on product quality. AOQ is a 
function of p. It has a maximum value which is called the average outgoing quality 
limit, AOQL. Tables are provided in MIL-STD 105D that give the factors with which 




















the AOQL values for specific sampling plans can be calculated. The AOQL is a value 
representing the worst case AOQ for a particular sampling plan. The value of p cor- 
responding to the AOQL is of particular importance. This subject will be addressed 
further in Chapter IV. 

A producer may depend on either his or a consumer's quality control department for 
post-production inspection of his products. This inspection results in nonconforming 
items being screened out and replaced by conforming items. Often, this is how sampling 
plans such as MIL-STD 105D are used in a manufacturing operation. However, this 
post-production type of inspection is wasteful because it allows time and materials to be 
invested in products which are not always saleable. 

One way to solve this problem, is to adopt an inspection strategy wherein a product 
can be inspected for quality in the early and middle stages of production. A second way 
of solving this problem is the adoption of a strategy of prevention. Clearly. it is much 
more effective to avoid waste by not producing unsaleable output to begin with. What 
follows in the next chapter is the basis for such a strategy. 
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HI. A DIFFERENT APPROACH - DESCRIPTION AND ANALYSIS 





In this chapter, the differences between statistical quality control and statistical 
process control are addressed. These differences lead to another outlook : wards quality 
assurance and an alternate method of dealing with it which is set forth in this chapter. 
The fundamental idea of the alternate method is the use of control charts by the pro- 
ducer to gain statistical control of the manufacturing process with the result ultimately 
being better quality in the product. Since a control chart is a form of contro! based upon 
sampling rather than 100 percent inspection, the quality of process output will depend 
not only on the control chart, but upon the characteristics of the process which cause 
its likelihood to produce defects to varv over ume. Thus, a model will be proposed to 
estimate the performance of the manufacturing process and a Ma....ov chain will be used 
to develop a measure of process output. Both of these will be presented in detail. 

A somewhat similar study was conducted bv Brugger and is detailed in his technical 
paper [Ref. 5]. Brugger did use a Markov chain approach; however, his study was fo- 
cused on modeling an acceptance sampling procedure vice a control chart procedure. 


A. STATISTICAL QUALITY CONTROL VERSUS STATISTICAL PROCESS 
CONTROL 

Quality assurance will now be approached in a different manner than in acceptance 
sampling. It will be addressed from the manufacturer's point of view. Aside from rec- 
tifving inspection, quality cannot be inspected into an item after it has been manufac- 
tured. Because of this we will focus our attention on the manufacturing process. If this 
different viewpoint seems to the reader to be a minor point, then it should be pointed 
out that it changes entir “* how an analyst will approach the problem and what tech- 
niques he must use to solve it. Instead of a statistical quality contro] (SQC) problem, 
we are now faced with a statistical process control (SPC) problem. SQC involves a large 
amount of inspection. This inspection is usually conducted on final product. However, 
at this point in the process it is too late to make corrections or adjustments to the 
process which can correct for the mistakes already made: although, it will certainly be a 
necessarv task if quality is to be improved in future items. SPC methods are bused upon 
the examination of finished ad semi-finished product at an early stage with some 
means of rapid and elfective feedback. Rapid feedback gives tighter control, saves add- 
ing value to defective items, saves time, and reduces the impact of defective material on 
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scheduling and resulting output. Effective feedback can only be achieved by the use of 
statistically based process quality control methods {Ref. 6]. SPC methods will allow the 
quality of the product to be controlled by allowing the process to constantly monitor 
itself. This means that the operator or technician in charge of the process examines the 
product as it comes through the process to make sure that it meets the required specifi- 
cations. He does this by using a sampling procedure which tells him exactly how often 
to take a sample and how many items to sample. The sampling must of course be con- 
ducted randomly. 

SPC methods are often verv simple when implemented in a manufacturing process. 
AS an example, suppose a quality problem arises within a process. The problem will be 
noticed reasonably quickly if the process is monitored correctly. After becoming aware 
of the problem, with minimum delay, the cause is determined and the problem is cor- 
rected. The procedure just described is a routine one and occurs frequently in a manu- 
facturing process. The use of SPC methods allows the frequent occurrence of procedures 
such as the one mentioned above to be dealt with promptly and efficiently. 

Problems with quality in a manufacturing process in which SQC procedures are be- 
ing employed are not always as straight forward. In SQC, a quality problem is some- 
times discovered after the product is out of the manufacturer’s factory and has been 
delivered to the customer. If through acceptance sampling procedures, the lot is re- 
jected, then not only has the manufacturer lost money but he has a defective product 
and it is conceivable that he may not even know why. 

The companies that have used SPC procedures over the years have found that 
qualivy costs related to SPC are usually Known and low. Once SPC procedures are im- 
plemented by a company, they tend to remain in use because they are found to be of 
considerable benefit. Unfortunately, the number of companies that are actually making 
Widespread use of SPC methods in their operation is relatively low when one considers 
the benefits that SPC methods offer. The major reason found for this low usage is a lack 
of knowledge, particularly among senior managers. Although these senior managers 
sometimes recognize quality as being an important part of corporate strategy, they do 
not appear .o know what effective steps to take in order to carry out the strategy. Too 
often quality is seen as an abstract property and uot as a measurable and controllable 


parameter. [Ref. 6] 
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B. CONTROL CHARTS AS A MEANS OF PROCESS CONTROL 

Assuming We do not want to conduct 100 percent inspection on the process we need 
to devise a sampling procedure or some other method that will tell us when the quality 
in our manufacturing process begins to deteriorate. In addition, we would like the 
method to be more versatile than just being able to tell us when we are producing ex- 
cessive defective products. We would like the method to be useful in helping us control 
the process, and a control chart is a widely used statistical method that is well-suited for 
this task. 

The control chart is considered to be a very important tool in statistical process 
control. Manufacturers readily accept the fact that the measured quality of their prod- 
uct is always subject to a certain amount of random variation which can be attributed 
to chance. Typically, some inherent random error is present in any process. Variation 
Within the process will, therefore. be inevitable. The control chart is a tool which gives 
the producer the capability to detect any excessive variation in the process which mav 
be attributed to assignable causes. Such assignable causes may be: differences among 
machines, materials, or workers, or a number of other factors related to the performance 
of the process, The reason for the excessive Variation is then determined and corrected. 
Thus, use of the control gives the producer the capability to diagnose and correct many 
production problems and often brings substantial improvements in product quality. 
Furthermore, by identifVing certain quality variations as chance variations, the control 
chart tells us When to leave a process alone, thereby preventing unnecessary aid frequent 
adjustments that tend to increase the variability of the process rather than decrease it. 

One versatile and widely used control chart is the p-chart. It is the control chart 
for proportion defective, or fraction defective, which is designated by p. A p-chart 
works as fellows. A sample of a predeternined size is taken at a predetermincd interval. 
The number of items in the sample that are defective, x . is then compared with the 
sample size, s . to form our ratio < , or our sample fraction defective which we will 
refer to asp. This value is then compared to the control limits on the control chart to 
get an idea of how the process is performing at the time of the sample being taken. The 
n-chart can be constructed using data that are already available for other purposes or 
can readily be made available. 

We will use the p-chart in our quality method. Before we can use the p-chart, we 
must have an estimate of the fraction defective of the process. We will assume that we 
are able to obtain this estimate without a great deal of difficulty. Now we have nearly 


all of the information we need to construct the p-churt. Simple statistical calculations 
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can be used to provide control limits that tell whether assignable causes for variation 
appear to be present or whether the variations from lot to lot are explainable by means 
of chance occurrence. We will concern ourselves only with an upper control limit, since 
a lower control limit on fraction defective is a self-defeating purpose. The upper control 
limit can be calculated as follows: 


UCL; = p +30; — [Ref 3]. 
We should be able to yet an estimate for p from past data on the process. We will 
denote this estimate as p. We will use the value for p as the centerline on the control 


chart. Also, since, p is based on the binomial distribution, we are able to estimate 


__ [pa=s) 
o = $ 


As can be seen a; is dependent on the sample size, s , which will vary from control chart 


o; as follows: 


to control chart. Thus, o; will also vary from control chart to control chart. 

For illustration purposes, suppose our sampling frequency, m , or the number of 
items produced between samples being taken, is 1000. Also, assume our sample size, 
s ,is 80 and our estimate p of the process fraction defective p is .015. We can now 


calculate o; as follows: 


ot 


__ [.015(.985) 
7 80 


0136. 


tt 


Now that we have a value for o; , we can calculate UCL; , or the upper control limit 


for the process as follows: 


UCL, = .015 + 3(.0136) 


0558 
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We are now able to proceed with the implementation of the control chart for our 
process. Figure 2 shows a control chart with the parameters defined as follows: 


oS = 0136, 


and 


mc 1000. 
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Figure 2, An Example of a p-chart 


One should keep in mind that the value for o; is heavily dependent upon the sample 
size, s . What the control chart is telling us in this example is that the process is out 
of control because it has one point above the upper control limit. As we discussed ear. 
lier, what should happen in this instance is that the source of the problem is determined 


as quickly as possible and then promptly corrected. 
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C. MEASURING THE EFFECTIVENESS OF CONTROL CHARTS 

Taking the concept of control charts one step further, we would like the consumer 
to be able to determine the effect that using a control chart at the end of a process will 
have on the quality of the output. The consumer can then judge for himself whether 
or not acceptance sampling is needed. Once again, the consumer’s primary concern is 
whether or not either procedure will function adequately to furnish product of accepta- 
ble quality. 

One commonly used measure of control chart performance is the chart’s operating 
characteristic curve or OC curve. A control chart’s OC curve gives the probability that 
(for a given fraction defective) the process is declared under control based on a sample 
result. Being declared under control is essentially the acceptance of the null hypothesis 
that the process fraction defective is p , the value based upon past data of the process. 
This probability as a function of the current process fraction defective p is denoted 
Pulp) . 

Another well-known measure of control chart performance within quality circles is 
the average run length curve or ARL curve. The ARL curve gives, for a particular 
fraction defecuve, the average number of samples that would have to be taken before the 
process was determined to be out of control. [Ref 1] Figure 3 shows an cxample of a 
p-chart OC curve and Figure 4 shows an example of an ARL curve. 
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Figure 3. p-chart OC Curve 
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Figure 4. ARL Curve 


Both the OC curve and the ARL curve are shown as functions of p , the process 
fraction defective. Furthermore, there is another less obvious relationship between an 
ARL curve and an OC curve. It can be shown that the ARL value is the mcan of a 


geometric distribution, and the average run Jength can be calculated as follows: 


= 1 
ARLQ) "T= pap » [Refi 1] 
Still, the consumer likely will not be satisfied with either of these methods since 
neither one will tell him the quality of received goods. The consumer would like to be 
able to evaluate the producer’s method of quality assurance to determine the quality 


level of the product he is receiving. With respect to the control chart procedure, the 

















question now becomes: how does the consumer judge a manufacturer’s control chart? 
Upon initial examination of a p-chart one notices the following: 

1, The sample size, 

2. The sampling interval, 

3. The upper control limit, which is a function of the centerline and sample size, and 


4. The steps to be taken by the producer in the event the process is determined to be 
out of control , i.e, when the sample proportion is found to be outside the control 
limits on the control chart. 
The customer’s concern is that none of above listed factors gives him any direct infor- 
mation about the quality of product that he will receive. 

In addition to these factors, it is clear that product quality leaving a final control 
chart will depend upon how the characteristics which are the primary cause for defects 
in the process change over time. Clearly, these characteristics will vary from process to 
process, and in some instances may be unknown or difficult to estimate. In the next 
section, we will give one simple model of a manufacturing process in which p varies 
over time. The model will be used to demonstrate the effectiveness of the control chart 
procedure, and to provide a basis for computing a measure of the quality of goods re- 
ceived by the consumer. 


D. A MODEL OF A MANUFACTURING PROCESS 

‘the next step in this approach to our problem is to devise a way of approximating 
or estimating the performance of the manufacturing process with respect to quality. The 
simple process we will model is one where a final control chart is emp] sed for overall 
fraction defective. We shall assume that when the control chart shows that the process 
is out of control, corrections are made before additional items are produced, and the 
process fraction defective returns immediately to a lower base value. This base value is 
the :esult of inherent, random errors within the process. 

When devising this model, one must include the effect that various factors have on 
the quality of a product. One of these factors is the aging of elements uscd in the process 
(machinery, tooling, etc.). Another factor is the random way in which defective products 
occur and how we record them; e.g.. by attributes or by variables. In our model, items 
are scored by attributes only; that is each item is classified into one of two categories: 
defective or non-defective. 

Another modeling problem that must be addressed is: what happens when the 
process is determined to be out of control or producing an excessive amount of defective 











products. Typically, some form of corrective action is taken to get the process back 
under control. At times, this corrective action may involve such measures as shutting 
the process down until the problem can be determined and corrected. Drastic measures 
such as this are usually not required since the probem can often be corrected without 
noticeably slowing down production. For modeling purposes, we will assume that when 
the process is back in control as a result of the necessary corrections being made, the 
process fraction defective immediately reverts back to a value which is a result of the 
occurence of inherent, random errors in the process. This seems to be a reasonable as- 
sumption when considering the reasons listed above for occurrence of defects in a 
product. The assumption of the process fraction defective reverting back to its inherent 
value or state along with the general nature of a manufacturing process allows us to 
model the process performance as a stochastic process. First, let us suppose we divide 
up the daily production output quantity of the process into # increments of size m . 
This can be done on a shift basis instead of a daily basis if the producer so desires. This 
gives us our sampling frequency. In our example we will arbitrarily choose 6 as the 
value of ma . Thus, we have six production periods per day or per shift and a sample 
Will be taken after every mi units of production. Let a be defined as the probability 
that an item produced by the process when under control is defective duc only to the 
inherent random error within the process. Whenever we have determined the process to 
be out of control. the model allows the process to become under contro! again by per- 
mitting the probability of producing a defective item to immediately revert back to a . 
The assumption is made that the problem causing the process to be out of control can 
be determined and corrected. Also, let us define 5 a» the rate per item produced at 
which the probability of occurrence of a defective item increases with respect to pro- 
duction. We now have the basis for our process model which is shown in graphical form 
in Figure 5, 

It should be pointed out that in a real world manufacturing process. the parameters 
aand 6 will not only be unknown but they will probably be verv diflicult to estimate. 
If they were known, then we would have no reason to employ a control chart on the 
process since we would know when to make the necessary adjustments to the process to 
keep it producing goods of an acceptable quality. Likewise. we would have no use for 
a sampling plan if we had available to us known values for a and b. By using ovi 


model we hope to hase a rough means of portraving a change in p over ume. 
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Figure 5. Graphical Representation of Model 


In Figure 5 we show the defect rate increasing linearly with the number of items 
produced. Sampling inspection will take place after every m= items are produced, 
When the process is determined to be out of control via use of the control chart, the 
defect rate immediately returns to the point designated in the figure by the coordinates 
(0,a). The process essentially starts over at this point. Moving along with the model, 
let Pa, equal the probability of concluding that the process is under control at level 
a+ibm where i is the period number and a+ibm_ is the process defect rate at the 
end of period i. The probability Pa, can be obtained from from the OC curve for the 
control chart. We can think of im as the number of items produccd since the last time 
the process was determined to be out of control and at which time the probability of 


producing a defective item was returned to (0, a). If we let P, equal the fraction 
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defective for the quantity of goods produced during period / , then we can calculate this 
value for each period as 


P,= a+ ae bm. 

As we noted earlier, we have built-in to our model the control chart decision rules. 
An important point is that without these decision rules in our model the process fraction 
defective, p , will increase linearly without bound. This is duc to the fact that the con- 
trol chart is the only measure we are using to decide whether or not the process is under 
control or not. However, as we indicated earlier the purpose of the model is simply to 
give us a rough idea of how a manufacturing process performs when a control chart is 
imposed on it. 


E. CALCULATION OF STEADY-STATE PROBABILITIES USING MARKOV 
CHAIN PROPERTIES 
Generally, if and when the manufacturing process reaches a level with defective rate 

a+6bm in our example, the probability that we conclude the process is under control 
will be extremely small. Therefore, we will assume that Pa, is equal to zero. This as- 
sumption gives us the capability to approximate the behavior of the process with a 
Markov chain with a finite number of states (six). A Markov chain can describe a 
process which moves from one state to the next in such a way that the probability of 
going to one of the next states only depends on the fact that it is presently in a given 
state and does not depend on anv prior states [Ref. 7]. We can summarize what we 
know thus far in Table 1. which follows. 











Table 1. SUMMARY OF MODEL PROBABILITIES 
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The six states in the model correspond to the ends of the six periods listed in Table 










Period Probability of a 
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1. State 1 is the end of Period 1 where the defect rate is a+bm. The Markov chain 
reaches State | as a result of the process being determined to be out of control. When 
this occurs, corrective action is taken on the process which takes the process defect rate 
back to a , and then with certainty to a@-+ bm at the end of Period !. 

Referring to Figure 5S, one can see the pattern which is a result of the model being 
a Markov chain. In State 1, the process can only transit to State 2 if it is concluded that 
the process is under control or to a if it is concluded that the process is out of control. 
Similarily, in State 2, the process can only transit to State 3 or to a. Let the matrix 
T , which follows, represent the transition probabilities. This 6 x 6 transition proba- 
bility matrix shows the probabilities of transitioning to states represented by the col- 
umns, given one is in a state represented by one of the six rows. Note that the transition 
probabilities are either 0, 1.0. or are from the OC curve from the control chart of the 


process. 














For an irreducible ergodic Markov chain, z,. or the long-run proportion of time that 


the process is in state j', can be calculated as follows: [Ref 7] 


and 


Using the transition probabilities, we can develop equations that can be used to solve 


for z, as follows: 


m = m(1 — Pay) + 2(1 — Pay) + 73(1 — Paz) + 23(1 — Pay) + 25(1 — Pas) + 7%, (3.1) 
R = mPa, ry (3.2) 
tO 2Pay ’ (3.3) 
ty = 73Pa; ’ (3.4) 
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ts = T4Paz, (3.5) 


% = m5Pas (3.6) 
and 

6 

Yay = 1.0. (3.7) 

jal 


We now have seven equations and six unknown variables, (we can take the Pa, values 
from the OC curve). We will work with equations (3.2) thru (3.7). If we add these 
equations together and perform some algebraic manipulation. we end up with the fol- 
lowing results: 


l—7, = 2,Pa, + 2,Pa,Pa, + 2,Pa,Pa,Pa, + 2,Pa,Pa,Pa,Pa, + 7,Pa,Pa,Pa,PaPa, , 


Which simplifies to 


] 


juan ee Pa, + PayPay + Pa, Pa,Pa, + Pa,Pa,Pa,Pay + Pa, Pa,Pa;Pa;Pas 


. (3.8) 


Equations (3.2) through (3.6) show that each xz, (i = 2,...,6) is proportional to 


z 


t,. We can develop the following equations for z, : 


Pa, 

ee eee ere 
: 1 + Pa, + Pa,Pa, + PaPa,Pa, + Pa,Pa,Pa,Pa, + Pa,Pa,Pa,Pa,Pa;' 
Pa Pa, 
QV une | 
: 1+ Pa, + PaPa, + PaPa,Pa, + Pa,Pa,Pa;Pa, + Pa,Pa-Pa,Pa,Pa,’ 
Pa, Pa,Pa, 
i a 
* 1+ Pa, + Pa,Pa, + Pa,PaPa, + Pa,Pa,Pa;Pa, + Pa,Pa,Pu;PaPa,’ 
Pa,Pa,Pa;Pa, 
Ts 


amen ar’ Pa, + PaPa, + PaPa,Pa, + Pa,Pa,Pa;Pa, + Pa,Pa;Pa,PaP ae 
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and 


Pa,Pa,Pa,Pa,Pa 


“= "T+ Pa, + Pa,Pa, + Pa,Pa,Pa, + Pa,Pa,Pa,Pa, + Pa,Pa,Pa,Pa,Pa, * 


Now that we have determined the long run proportion of time the process spends 


in each state, and Table | shows the fraction defective in each state, we can calculate the 
fraction defective of the process as determined by the model. We will denote the model's 
resulting value for the process fraction defective as p’. 

We can compute p’ as follows: 


p’ = nlat + bin) + a(a+ 4 bin) + 2(a+ 7 bm) + rat = bi) 


+ rat? bm) + adatt dm). 
If we collect the common terms, we end up with the following equation for fraction de- 
e fective as computed ‘vith the use of our model: [Ref. 7] 


6 


pa ad bm >x( 45+). (3.9) 


f=1 


The expression for p’ in Equation (3.9) provides a wav to compute a value for the 
output fraction defective of the simple process we have modeled. The six 2; values 
may be obtained from solution of Equation (3.8) for 2, , and successive solution of 
Equations (3.2) through (3.6) for 2, through , . Solving these equations requires five 
values: Pa,, Pa, Pa,, Pa, and Pa, . These values can be obtained from the control chart 

- for the manufacturing process. 

An important point to note is that p’ may not accurately reflect the level of quality 

in the product as it ends up in the consumer’s warehouse. This is because the model 


does not consider how the quality level in the product is affected by shipping . onloading, 











oMoading, etc.. Despite the model’s limitations, it gives us a rough estimate of a proc- 
ess’ fraction defective when a control chart is employed. 

The end result of the model is the value of the manufacturing process’ output frac- 
tion defective, or more simply put, an estimate of the level of quality in the product 
produced by the process when there is a control chart at the end. This value p’ can 
be thought of as a means of measuring the effectiveness of the control chart procedure. 
Furthermore, this value can be compared to the AOQ value of an acceptance sampling 
procedure to get an idea of which method is more effective. This will be addressed in the 
next chapter. 


F. ANALYZING SOME RESULTS OF AN IMPLEMENTED CONTROL CHART 
AT THE END OF A PROCESS 

It should be clear that the parameters of the model can be changed as necessary in 
order to estimate the performance of a variety of manufacturing processes. The program 
listed in Appendix A was written so we would have a means of testing the model with 
respect to its ability to estimate the output of the process with a contro] chart. The 
program allows the user to judge the sensitivity of the model with respect to variation 
of its parameters a@ and 6, as well as the contro] chart characteristics UCL;, sample 
size s, and sampling interval mm. By varying these factors we are able to get a better 
understanding as to which ones are the most critical in computing the process output 
average fraction defective. p’. 

The program was applied to cases where a varied from 0 to 0.0075, b from 
0.900001 to 0.000005, p from 0.001 to 0.01, sample size s from 50 to 150, and sampling 
frequency mm from 500 to 2000, 

Some results of the computer program (see Appendix B), suggest that. over the 
range of parameter values used, the estimate of average fiaction defective fluctuated the 
most as the parameters 6 and m were varied. Hence, & and m or their product appear 
to have the most influence on the estimate of the output quality of the process, This 
seems logical when one recalls that mm is the sample frequency and bis the rate at which 
the process performance deteriorates. Ifa process has a high rate of deterioration, then 
it is clear that its average fraction defective will be higher as a result even with the use 
of a control chart. Still, it will be significantly lower when using a contro! chart than 
when not using one. Intuitively, this will be true of any manufacturing process in a real 
world environment and, as explained earlier, it is true by design in our model. This is 


because our model cannot operate without a control chart. Also. if quahty problems 

















are occurring within a process, the more frequently samples are taken, the greater the 
frequency of detecting and correcting those problems becomes. 

Taking our analysis one step further, we will now focus our attention on a portion 
of the data. For the cases where the sampling frequency is m = 1000 ( one third of the 
data), using the data in Appendix B, we further find that the process fraction defective 
is most dependent on the process characteristics, a and b , and the sampling frequency 
m.. This was determined by obtaining a linear fit for the data. In particular, the process 
fraction defective after one period. (a + bm) reflects both the process characteristics and 
the sampling frequency. Using this as one variable and the control chart rule (UCL) as 
the other, we applied a least squares fit to our data. The result for the $1 data points 
Where m= 1000 was 


p’ = 0.0005 + 1.0094(a + bm) + 0.0733(UCL) . 


This reflects, for this set of data, the clear influence of process characteristics and sam- 
pling frequency. Some values of interest for the m= 1000 set of data are shown in 
Table 2 which follows. Note that this is not a statistical regression analysis, but simply 
an effort to fit a linear function to a set of data. 


Table 2. RESULTS OF LEAST SQUARES FIT TO DATA 


a 
Value Value 
0133 


Resid 





Maximum 









As stated earlier, the purpose of the computer program shown in Appendix A was 
simply to demonstrate the effectiveness of a final control chart on process output as a 
function of process parameters and contro] chart characteristics. It is believed that if 
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one examines the output, its effectiveness will become apparent. After close examina- 
tion, one can also get an understanding of the sensitivities of the model to specific pa- 
rameters. This information can be useful to a quality manager when he is attempting 
to improve his operation. 

Although, the computer program provides a means of evaluating the effectiveness 
of the control chart procedure, what is really desired is a method of comparing, from the 
consumer's point of view, reliance on a final control chart to acceptance sampling. This 
issue is addressed in Chapter IV. : 
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IV. COMPARISON OF THE RELIANCE ON ACCEPTANCE SAMPLING 
VERSUS A CONTROL CHART PROCEDURE 


In this chapter, we will look at various means of comparing acceptance sampling 
With statistical process control, wherein acceptance sampling would not be employed. 
First, a qualitative comparison of the corresponding costs of each method will be per- 
formed followed by a comparison of the effectiveness of each method. 


A. COST COMPARISON 

The cost comparison will be approached from the perspective of the base 
commander. From his viewpoint, anv alternative that reduces jis costs and maintains 
or increases current levels of effectiveness will attract his interest. Since the proposed 
alternative may reduce the the number of quality assurance personnel required to be 
staffed at the base, the base commander will at least reap the amount of savings resulting 
from this. However, the biggest savings in using a method of quality assurance involving 
control charts vice acceptance sampling is what results from the elimination of sampling 
inspection and rectifving inspection. 

An outsider may argue that even though in this case the Naval Weapons Support 
Center is saving moncy, the overall cost to the Department of Defense is the same. 
While it is true that there is a shift in cost from the naval base to the defense contractor, 
there still may be a significant savings involved with the control chart procedure. First 
of all, the defense contractor should have no need to hire additional] personnel when he 
implements this method. This is because it is necessary and fundamental to the success 
of the process control method that the technicians and operators on the production line 
should be the personnel who are carrving out the methods of the control chart proce- 
dure. This is where training of personnel is important. It may be necessary to train all 
of the technicians and machine operators associated with the manufacturing process in 
the fundamentals of statistical process control and the use of contro! charts. The major 
impact of training costs will be felt in the beginning while the changeover to rcliance on 
the control chart procedure is being made. Training will still have to be conducted, 
though to a lesser degree. on a continual basis. This will be one of the two primary costs 
of changing the emphasis from an acceptance sampling procedure to a control chart 
procedure. The other cost is associated with instances when it is necessary to stop the 
process to determine the reason for the process being out of control. This resulting 
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downtime in production can be verv costly and can quickly cut into the defense con- 
tractor’s profit margin. 

The cost considerations discussed above are difficult to quantify. Why these costs 
which primarily affect the manufacturer are of concern to the base commander is that 
he will see them reflected in his contract with the manufacturer. Due to the factors 
mentioned above, the base commander may sce a reduction in his costs when the con- 


tract is written if a control chart procedure is used. 


B, EFFECTIVENESS COMPARISON 

Having concluded that elimination of acceptance sampling may result in a cost 
savings, We now to focus our attention on the effectiveness of a control chart procedure 
as opposed to an acceptance sampling procedure. One wav of examining the effective- 
ness of each procedure is to look at the effect of each on the AOQ or average outgoing 
quality of the process and the quality components. 

With the concept of AOQ, we will assume a scheme of 100 percent inspection of 
rejected lots and that any defective items discovered within the rejected lot will be re- 
placed by good items. Another assumption is that the 100 percent inspection of the re- 
jected lots will ultimately result in 100 percent quality of those lots. Finally. in our 
analysis, we will assume that the number of defects in our sample is binomially distrib- 
uted with a mean equal to sp. or the product of the sample size s and lot quality p. 
We will now show how average outgoing quality is related to p. 

Recall that in Chapter I, the AOQ for a process in which a sampling plan and 100 


percent rectification are implemented was defined as follows: 


Pa(p(L — s)) 


AOO = 7 


For illustration purposes, we will use the same MIL-STD 105D sampling plan that 


we used earlier as an example. Its parameters are as follows: 


lot size (L = 100), 
sample size (s = 80), 
acce, table quality level (AQL = .015) and 


acceptance number (c = 3). 
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(Note: MIL-STD 105D uses the following letter designations: N 

for lot size (vice L) and n for sample size (vice s). Different 
designations are used in this thesis to clarify the difference in the 
parameters used in the model and the parameters used in MIL-105D.) 


Using various values for the quality of incoming product, p, we summarize the re- 
sulting AOQ values for the different values of p in Table 3. The values for Pa were 
taken from the OC curve for this plan shown earlier in Chapter I]. 


Table 3. AOQ VALUES FOR A MIL-STD 105D SAMPLING PLAN WITH 
AQL=.015 
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As can be seen the maximum AQQ value is about .0223, This AOQ value is known 
as the AOQL. or Average Outgoing Quality Limit. The corresponding value of p at 
which this maximum AOQ occurs is .035. Simply put, the process with use of this par- 
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ticular sampling plan is ultimately producing its worst level of quality, or 
AOQL = .0223, when p, the lot quality, is equal to .035. In effect, we can consider 
this to be the worst case performance of this MIL-STD 105D sampling plan when rec- 
tifying inspection is employed. 

Under the proposed method of using control charts to gain statistical control of the 
process, our argument is a very simple one. The process long-run average performance 
p Will simply be equal to the estimate p’ of the process average fraction defective ob- 
tained from the modcl. Figure 6 shows the effect of each method on the AOQ of the 


process. 
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Figure 6. Comparison of the Effect of the Two Methods on AOQ of Process 


It can he seen from the figure that at times, MIL-STD 105D performs better and 


at other times it does worse. In this example the worst case performance is better with 

















the control chart method. The MIL-STD 105D AOQ curve is for the same sampling 
pian we have been using as an example throughout this thesis. Once again its parame- 
ters are shown below: 

L= 1000, 

s= 80, 
and 


AQL= 015. 


The two dashed lines in the figure representing the upper and lower bounds on p’ for 
the proposed method were computed using the following parameters in the model: 


n= 6, 
ne= 1000, 
a= 0075, 


b= .Q00001 (for the lower bound). 


b= .000005 (for the upper bound). 


s= 80, 


and 


in 


p= 01 


As we have already established, p’ is dependent upon the parameters uscd in the 
model. However, in this example the values for the model’s parameters were chosen so 
that the two methods will be operating on processes with roughly the same character- 


istics. The same value that was used for lot size in the military standard was chosen as 
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sampling frequency in the model. Also, it was decided it was best to use the same sample 
size in each method. The inherent probability that the process will produce a defective 
item a was arbitrarily chosen to be half of the AQL value. It is believed that for most 
manufacturing processes this is an overstatement of the tendency for the process to 
produce a defective item. The rate at which this probability of producing a defective 
item increases with time 6 was chosen to take on two values for illustrations purposes 
due to its uncertainty. The lower bound value represents what is believed to be a rea- 
sonable deterioration rate for the process; whereas, the upper bound represents an ap- 
proximation of a worst case scenario. It is difficult to compare the two methods on a 
equal basis; however, it is believed that the estimates made here are conservative and 
provide for a rough, but reasonable comparison of the two methods. 

Hence, we conclude that by using this means of comparison, it is reasonably evident 
that implementation of the proposed method of reliance on a control chart procedure 
should be given consideration because it will likely result in a cost savings. Also, our 
analysis has shown that it some cases it can be at least as effective as the current method 
of acceptance sampling which is being used. 

In the following chapter. we will summarize our study and make further conclusions 
with regard to what we believe the effect of a control chart procedure similar to the one 
presented here will have on a manufacturing process if implemented. 




















V. SUMMARY AND SUGGESTIONS FOR FURTHER STUDY 


The purpose of this thesis was to conduct research into the area of quality assurance 
and to determine whether or not there might exist a feasible alternative to ensuring 
quality in military systems as compared to the methods being used today. Specifically, 
the Commanding Officer of the Naval Weapons Support Center in Crane, Indiana, 
would like to know if it would be feasible for him to rely only on the manufacturer’s 
operation and quality assurance methods to ensure an acceptable level of quality in the 
product he is receiving from the manufacturer, and how he could evaluate these meth- 
ods. Since there seems to be a slow but sure shift in the manufacturing industry towards 
Statistical process control, it was decided that this thesis would concentrate on an alter- 
native of this nature. 

Acceptance sampling conducted by the consumer requires that the military agency 
staff a group of personnel knowledgeable in quality assurance methods. The method 
presented in this thesis, consistent with the fundamental concept of statistical process 
control, shifts the emphasis of product quality back to the manufacturer where ideally 
problems in quality are detected early and corrected on the spot. In order for the cus- 
tomer to shift from the reliance on acceptance sampling to reliance on control charts at 
the manufacturing operation, is necessary to show that upon consideration of the pro- 
ducer’s quality program, the product he will receive without acceptance sampling is sat- 
isfactory. To do this, an effective measure is needed to compare the alternatives of 
acceptance sampling and control charts. 

As a way of measuring the effectiveness of the control chart on the output of a 
manufacturing operation, a model of a manufacturing process was developed. The pri- 
mary assumption made in constructing the model was that a linear deterioration rate 
over time existed with respect to quality. Another assumption made was that when the 
process is determined to be out of control, the problem can be determined and corrected. 
Upon correcting the problem, the value for the rate at which defects occur returns im- 
mediately to a value that is due only to the random error present within the process. 
The model requires inputs for sample frequency; sample size; the probability of a defect 
occurring due to inherent or random error present within the process; the rate at which 
the probability of the occurance of a defect increases with time; and the centerline value 
for the control chart to be used. The output p’ of the model is an estimate of the 
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process average output fraction defective when the control chart is in use. The 
paramters in the model that were determined to be the most influential in this estimate 
were the sample frequency and the rate at which the probability of an occurence of a 
defective product increased with respect to time. 

Although, the main purpose of the model is to show one way to evaluate the effec- 
tiveness of a control chart, a quality manager may be able to use it as a decision aid if 
control charts are being employed in his operation and if he believes his manufacturing 
process possesses the necessary characteristics (i.e., linear deterioration rate, etc.). The 
manager will have to be able to estimate within reasonable accuracy the parameters in 
the model as they apply to his process in order for the model to provide him with useful 
information. The quality manager will be able to derive his estimates from existing data 
concerning past performance of the process. Even if useful data does exist on the per- 
formance of the process, the parameters in question will probably be difficult to estimate 
due to the uncertainty of their behavior, as will probably be borne out by the data. 
Clearly, if the manager were to have at his disposal estimates of these parameters which 
involved little or no uncertainty then he would have no use for a control chart to begin 
with. Nonetheless, if the quality manager is able to obtain estimates of a and b then 
he may be able to use the model as a decision aid to assist him in improving the effi- 
ciency of his process. 

In Chapter IV, the argument was made that even if reliance on a control chart 
without acceptance sampling resulted in approximately the same level of quality in the 
product, it would result in a cost savings because it would reduce the need for staffing 
of quality assurance personnel by the customer. We also argued that the proposed 
procedure would reduce costs because a manufacturer is much better equipped to per- 
form a quality assurance operation than the consumer. For statistical process control 
to be effective, the quality assurance methods miusz be performed by the manufacturer. 

There is much room for improvement in manufacturing methodologies and quality 
assurance techniques. It is recommended that this research be taken further by ap- 
proaching it slightly differently. An alternative approach to the problem is to model the 
manufacturing process using two rates of average fraction defective and the time be- 
tween the process changing between the two would be exponential. There are many 
manufacturing systems that would lend themselves to this tvpe of model. 

It is sincerely hoped that the methods and procedures presented in this thesis will 
be useful in helping the Naval Weapons Support Center and other Department of De- 
fense agencies in dealing with quality assurance issues. 
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APPENDIX A. 


Program to Generate OC Curve Probabilities 


INTEGER M(1: 10) ,MM,N,S(1: 10) ,SS,A4A,BB,C, ICOUNT, PP,NLINE 
REAL PA(1:10),PI(1: 10) ,FRADEF(1: 10) ,SUM,A(1: 10) ,B( 1: 10) 
REAL PTOTAL,CL( 1: 10) ,UCL(1: 10,1: 10) ,DENOM(1: 10) ,DSUM 
REAL PPRINE(1: 10,1: 10,1: 10,1: 10,1: 10) 


ve ve de te de ve ve te te ode ve de ove ve ve ve ove ve te ke ove te ve ve te te ve ve te ee te ke oe 
GIN CCL LUEL TENE MTATE GE ELULE TELE RMUEL AE EG, . 
// The purpose of this program is to demonstrate the // * 
// effectiveness of a control chart procedure when it is // * 
// imposed upon a manufacturing process. This is done by // * 
// evaluating a model of the process. Given the necessary // * 
// inputs, the OC curve probabilities are calculated. // * 
// Since, the process is modeled via the use of a markov // * 
// chain, the steady-state probabilities are then // * 
// calculated, These calculations allow for the // * 
// computation of an estimate of a fraction defective // * 
// for the process in question. // * 
1 es 
{TI RAMU REA RARER REA EREA SEALS LALA REALE REAL RUAEAEAESEAEAEA : 


Prompt the user for data entry. 


PRINT *, 'To check the sensitivity of the model to each parameter’ 
PRINT - ‘enter the number of different values you would like to' 
PRINT om "enter for each parameter. 


PRINT 

READ *, NUM 

PRINT * 

PRINT *, ‘Enter the value for the number of periods, n.' 


PRINT *, ‘This value will not vary. To change it you must’ 
PRINT *, ‘run the program again.‘ 
PRINT a 
READ **, N 
PRINT * 
DO 5 I = 1,NUM 
PRINT * 
PRINT *, 'Enter the value for the sampling frequency, m.' 
PRINT * 
READ *, MCI) 
PRINT * 
PRINT *, ‘Enter the value for the sample size, s.' 
PRINT * 
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aaan 


aqaaaa 


5 


READ *, S(I) 
PRINT * 


PRINT *, ‘Enter the value for a, the probability the process' 


PRINT *, ‘will produce a defective item when it is performi 3' 
PRINT *, ‘at peak efficiency. ' 
PRINT * 

READ *, ACI 

PRINT # 

PRINT *, 'Enter the value for B which is the rate at which' 

PRINT *, 'the probabilty of producing a defective item’ 

PRINT *, ‘increases with respect to time. ' 

PRINT * 

READ *, BCI 

une 

PRINT *, 'Enter the value for the centerline on the control’ 

PRINT *, ‘chart. 

PRINT * 

READ *, CL(I) 

PRINT * 
CONTINUE 
ICOUNT = 0 
NLINE = 0 
DO 60 MM 


Calculate the upper control limit for the control chart 


UCL(PP,SS) = CL(PP) + 3 * SQRT(CL(PP) * (1 - CL(PP)) 
/ S(SS)) 
C = §(SS) * UCL(PP,SS) 


Use the subroutine POISS to calculate the OC curve probabilities. 


10 


DO 10 I = 1,N 
CALL POISS(A,B,1,M,S,PA,AA,BB,MM,SS,C) 
CONTINUE 


Calculate the steady-state probabilities. 


15 
20 


PTOTAL = 
DSUM = 1. 
DO 20 J = N,2,-1 

DENOM(J) = 1.0 

DO 15 I = J-1,1,-1 

DENOM(J) = DENOM(J) * PACT) 

CONTINUE 

DSUM = DSUM + DENOM(J) 
CONTINUE 
PI(1) = 1 / DSUM 
PTOTAL = PTOTAL + PI(1) 


0 
0 
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DO 25 I = 2,N 
PI(I) = PI(I-1) * PA(I-1) 
PTOTAL = PTOTAL + PI(1) 
25 CONTINUE 


Scale the probabilities to ensure that they sum to one. 


DO 30 I = 1,N 
PI(I) = PI(I) / PTOTAL 
30 CONTINUE 


Calculate the fraction defective for each period. 


SUM = 0 
DO 35 I = 1,N 
FRADEF(I) = PI(I) * ((2*I - 1) / 2) 
SUM = SUM + FRADEF(I) 
35 CONTINUE 


Calculate the average fraction defective for the process. 


PPRIME(PP,BB,AA,SS,MM) = ACAA) + (B(BB) * M(MM) * SUM) 
40 CONTINUE 
45 CONTINUE 
50 CONTINUE 
55 CONTINUE 
60 CONTINUE 


Print out a definition list of the parameters. 


WRITE (4, 65) 

65 FORMAT ('1',/,15X, Definition of parameters’ , 
+/, 1X,'m = sampling frequency',/,1X,'s = sample size',/,1X, 
tia = inherent probability of the process producing a defective’ ’ 
oe item',/,1X,'b = rate at which the probabilitiy of producing’, 
t+' a defective’ »/,9X,'item increases over time’ s/s 1X, 'CL =', 
+' desired centerline value on control chart',/, 1X, ‘UCL, =', 
+' computed upper control limit for control chart’ Js 1X,'p = 
+’ resulting estimate of the process average fraction defective’ 11) 


Print a summary of parameter values for which the program was run. 


WRITE (4, 70) | NUM ,N 
70 FORMAT (10x, ‘The number of different values entered for each',/, 
+10X, ' parameter (excluding n) was: ',11,/,10X,'The value entered’, 
+' for n was: ',12) 
DO 80 I = 1,NUM 
WRITE (4, 75) M(I),S(1I),ACI),BC1I),CLCI) 
75 FORMAT (2X,'The following parameter values were entered: | des 
5¥,'M = "14, 4X, S = ',13,4X,'A = ',F5.4,4x,'B = Z 
F7.6,4X, "CL = ' ,F5.4) 


++ 
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80 CONTINUE 


Printout column headings for table. 


WRITE (4,85 


85 FORMAT (2X ‘case’ ,SX,'m" ,7% 's' ,6X,'a' ,9X,'b! , 9X, 


+ 'CL' 7X, 'UCL',7X,'p'/, 2K, !---=! 2X, '=--9" 5X, 
+ fore y3X, "ences "3X, !ooceeenee "AK, beeen! 3X, ence wee ey >: 
- eeesens /) 
Print out all calculated values. 
DO 120 MM = 1,NUM 
DO 115 SS = 1,NUM 
DO 110 AA = 1,NUM 
DO 105 BB = 1,NUM 
DO 100 PP = 1,NUM 
ICOUNT = ICOUNT + 1 
NLINE = NLINE + 1 
IF (NLINE .LT. 40) GOTO 90 
WRITE (4,65) 
WRITE (4,85) 
NLINE = 0 
90 WRITE (4,95) ICOUNT,M(MM) ,S(SS) ,A(AA) ,B(BB) ,CL(PP), 
+ UCL( PP, SS) , PPRINE( PP, BB, AA, SS ,MM) 
95 FORMAT (2X,14,2X,15,4X,14,3X,F5.4,3X, 
+ F9.8,3X,F5.4,3X,F5.4,3X,F7. 6) 
100 CONTINUE 
105 CONTINUE 


110 CONTINUE 
115 CONTINUE 
120 CONTINUE 
STOP 
END 
SUBROUTINE POISS(A,B,1,M,S,PA,AA,BB,MM,SS,C) 
INTEGER K,N(1: 10) ,S(1: 10) ,MM,SS,AA,BB,C 
REAL*4 ERROR, TOTAL,EPSLON,NU(1: 10) ,PAC1: 10), AC 1: 10) ,B( 1: 10) 


en cate nae a ome creme nee rr ten te re SN fe ce i AN Se eye ha-mm me Sa S08 sat Pry euRa-enemeney seem ace ce nmiemn mann nsec =e 


Uh ORE TA ERLE LTT CRU MEOREERDULELSE TAPED TEE EPAGLIEUEUL GUE 


// This subprogram will compute poisson probabilities // 
// by making use of the following equation: // 
// // 
// P(X = K) = ((EXP(-MU)) * (MU**K))/K ! , // 
// // 
// Where MU = a parameter based on the sample size and // 
// process quality. // 


// // 
LTITTTTTTT TTT TT TT TT TTT TTT TT TTT TAT 


ee ea a ee eee et ee A Se a RE i erry i oon a Segue new went Ps nen NS me ene ee ne eee are Sa Or mn a mr oe TE 
RS en ee re ee re a Ee OT TS NT A TT A Se SA AT TT ee Ee RE ES aS 
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MU(I) = S(SS) * (A(AA) + (I-1) * BCBB) * MCMM)) 


C 
C 
C Initialize variables 
. C 
K=0 
Cc 
; C When K = 0, to avoid division by 0, the 'prob' for K = 0 is simply 
C EXP( -MU) 
C 
PA(I) = EXP( -MU(I)) 
TOTAL = 0. 
C 
C 
C Begin summing the probabilities to find the cumulative probability. 
C 


TOTAL = TOTAL + PACT) 
5 IF ((K .GE. C ) .OR. (TOTAL .GT. .99999)) GOTO 10 
K=K+1 
Calculate the probability of exactly K occurrences. 


PA(I)= (PACI) * MUCI)) / K 


agqgand anaand 


Continue summing the probabilities to find the cumulative probability 


. TOTAL = TOTAL + PA(I) 
GOTO 5 





Set Pa(i) equal to the cumulative probability and return to main 
program 


FPOAIMAIAQ 


0 PA(I) = TOTAL 
END 











APPENDIX B. 


Results of Computer Program 


DEFINITION OF PARAMETERS 


m = SAMPLING FREQUENCY 

s = SAMPLE SIZE 

a = INHERENT PROBABILITY OF THE PROCESS PRODUCING A DEFECTIVE ITEM 

b = RATE AT WHICH THE PROBABILITY OF PRODUCING A DEFECTIVE 
ITEN INCREASES OVER TIME 

cl = DESIRED CENTERLINE VALUE ON CONTROL CHART (ROUGH ESTIMATE OF 
THE PROCESS AVERAGE FRACTION DEFECTIVE) 

UCL = COMPUTED UPPER CONTROL LIMIT FOR CONTROL CHART 

p = RESULTING ESTIMATE OF THE PROCESS AVERAGE FRACTION DEFECTIVE 


THE NUMBER OF DIFFERENT VALUES ENTERED FOR EACH 


PARAMETER (EXCLUDING n) WAS: 
THE VALUE ENTERED FOR n WAS: 6 


THE FOLLOWING PARANETER VALUES WERE ENTERED: 
m= 500 s= 50 a = .0000 b = .000001 


THE FOLLOWING PARAMETER VALUES WERE ENTERED: 
m = 1000 s = 100 a = .0050 b = .000002 


THE FOLLOWING PARAMETER VALUES WERE ENTERED: 
m = 2000 s = 150 a = .0075 b = .000005 


Case m s a b cl UCL 
1 500 50 . 0000 . 000001 - 0010 0144 
2 500 50 . 0000 - 000001 . 0025 . 0237 
3 500 50 .0000 .000001 .0100 = ,0522 
4 500 50 =—,0000~=—s «000002 )=—. 0010) =. 0144 
5 500 50 . 0000 . 000002 0025 0237 
6 500 50 =6.0000 +=.000002 +~.0100~ .0522 
7 500 50 .0000 .000005 .0010 .0144 
8 500 50 .0000 ~~ .000005 , 0025 . 0237 
9 500 50 .0000 = .000005 £0100 38,0522 

10 500 50 )=6, 0050 =. 000001 =. 0010 0144 
li 500 50 .0050 .000001 #£.0025 0237 
12 500 50 =6©.0050 =©.00000i1 .0100 .0522 





cl = .0010 


el] = .0025 














500 
500 
500 
500 
500 
500 
500 
500 
500 
500 
500 
500 
500 
500 
500 
500 
500 
500 
500 
500 
500 
500 
500 
500 
500 
500 


. 0050 
. 0050 
. 0050 
. 0050 
. 0050 
. 0050 
0075 
. 0075 
.0075 
.0075 
0075 
. 0075 
.0075 
0075 
.0075 
. 0000 
- 0000 
. 0000 
- 0000 
. 0000 
- 0000 
- 0000 
. 0000 
. 0000 
. 0050 
. 0050 


. 000002 
. 000002 
. 000002 
- 000005 
. 000005 
- 900005 
. 000001 
. 000001 
- 000001 
. 900002 
- 000002 
- 000002 
. 000005 
. 000005 
- 900005 
- 000001 
. 000001 
. 000001 
. 000002 
. 000002 
. 000002 
- 000005 
. 000005 
. 000005 
. 000001 
. 000001 
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. 0010 
- 0025 
. 0100 
- 0010 
» 0025 
. 0100 
. 0010 
0025 
. 0100 
. 0010 
. 0025 
. 0100 
. 0010 
0025 
. 0100 
. 0010 
. 0025 
. 0100 
- 0010 
. 0025 
. 0100 
- 0010 
- 0025 
. 0100 
. 0010 
. 0025 


0144 
- 0237 
- 0522 
. 0144 
. 0237 
- 0522 
0144 
. 0237 
. 0522 
.0144 
. 0237 
» 0522 
0144 
. 0237 
- 0522 
. 0105 
»0175 
- 0398 
~ 0105 
.0175 
. 0398 
- 0105 
.0175 
. 0398 
. 0105 
.0175 





. 0066 
. 0074 
. 0075 


0084 


. 0106 
O11 
- 0082 
. 0086 
. 0087 
. 0088 
. 0098 
. 0100 
0104 
. 0128 
. 0135 
. 0012 
0012 
. 0012 
. 0024 
. 0024 
. 0025 
.0055 
. 0055 
» 0062 
. 0061 
. 0061 








DEFINITION OF PARAMETERS 


SAMPLING FREQUENCY 

SAMPLE SIZE 

INHERENT PROBABILITY OF THE PROCESS PRODUCING A DEFECTIVE ITFM 

RATE AT WHICH THE PROBABILITY OF PRODUCING A DEFECTIVE 

ITEM INCREASES OVER TIME 

cl = DESIRED CENTERLINE VALUE ON CONTROL CHART (ROUGH ESTIMATE OF 
THE PROCESS AVERAGE FRACTION DEFECTIVE) 

UCL = COMPUTED UPPER CONTROL LIMIT FOR CONTROL CHART 


p RESULTING ESTIMATE OF THE PROCESS AVERAGE FRACTION DEFECTIVE 
Case m s a b cl UCL p 
39 500 100 .0050 .000001 .0100 = .0398 0062 
40 500 100 ~=6.0050 +=.000002 .0010 .0105 .0070 
41 500 100 =©.0050 .000002 .0025 .0175 .0070 
42 500 100 .0050 .000002 .0100 .0398 .0075 
43 500 106. —-«.0050)=. 000005 =3=—.0010 Ss. 0105s «. 0093 
44 500 100 =.0050 .000005 .0025 .0175 .0093 
45 500 100 =.0050 .000005 .0100 .0398 .0111 
46 500 100 . 0075 . 900001 .0010 .0105 . 0084 
47 500 100 0075 000001 .0025  .0175 . 0084 


48 500 100 .0075 .00000i1 .0100 .0398 .0087 
. 000002 =.0010 .0105 .0092 
50 500 100 =.0075 .000002 .0025 .0175 .0092 
51 500 100 =©.0075 .000002 .0100 .0398 .0100 
52 500 100) =.0075 =. 000005 = «60010 «6 0105). 0112 
53 500 100 = .0075 =.000005 =.0025) 60175. 0112 
54 500 100 =.0075 .000005 .0100 .0398 .0134 


> 
w 
Lea) 
oOo 
Co 
rary 
oO 
oO 
oS 
oO 
~I 
uw 


55 500 150 .0000 .000001 .0010 .0087 .0012 
56 500 150 =. 0000=S . 000001) =.0025 =. 0147. 0012 
57 500 150 =.0000 .000002 .0100 .0344 .0012 
38 500 150 .0000 .000002 .0010 .0087 .0024 
59 500 150 =.0000 =.000002 = .0025 .0147 .0025 
60 500 150 .0000 .000002 .0100 .0344 .0025 
61 300 150 .0000 .000005 .0010 .0087 .0049 
62 500 150 = .0000 .000005 .0025 .0147 .0058 
63 500 150 .0000 .000005 .0100 .0344 .0062 
64 500 150 .0050 .000001 .0010 .0087 .0059 
65 500 150 .0050 .000001 .0025 .0147 .0062 
66 500 150 =.0050 = .000001 .0100 .0344 .0062 
67 500 150 .0050 .000002 .0010 .0087 .0066 
68 500 150 .0050 .000002 .0025 .0147 .0072 
69 500 150 .0050 .000002 .0100 .0344 .0075 
70 500 150 .0050 .000005 .0010 .0087 .0082 
71 500 iso). 0050=—.000005 Ss. 0025S «6 0147S. 0098 
72 500 150 .0050 .000005 .0100 .0344 .0112 
73 500 150 .0075 .000001 .0010 .0087 .0082 
74 500 150 =©.0075 = .000001 .0025 .0147 .0085 
75 500 150 »=.0075 =.000002 =.0100) =.0344 =. 0087 
76 500 150 .0075 .000002 .0010 .0087 #.0087 
77 500 150 .0075 .000002 .0025 .0147 .0095 
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DEFINITION OF PARAMETERS 


n = SAMPLING FREQUENCY 

s = SAMPLE SIZE 

a = INHERENT PROBABILITY OF THE PROCESS PRODUCING A DEFECTIVE ITEM 

b = RATE AT WHICH THE PROBABILITY OF PRODUCING A DEFECTIVE 
ITEM INCREASES OVER TIME 

: cl = DESIRED CENTERLINE VALUE ON CONTROL CHART (ROUGH ESTIMATE OF 
THE PROCESS AVERAGE FRACTION DEFECTIVE) 

UCL = CONPUTED UPPER CONTROL LIMIT FOR CONTROL CHART 





P RESULTING ESTIMATE OF THE PROCESS AVERAGE FRACTION DEFECTIVE 
Case m Ss a b cl UCL P 
78 500 150 .0075 .000002 .0100 .0344 «.0100 
79 500 150 .0075 .000005 .0010 .0087 # .0100 
80 500 150 .0075 .000005 .0025 .0147 .0116 
81 500 150 .0075 .000005 .0100 .0344 .0136 
82 1000 50 .0000 .000003 .0010 .0144 .0022 
83 1000 50 .00C0O0 .000001 .0025 .0237 .0025 
84 1000 50 .0000 .000001 .0100 .0522 .0025 
85 1000 50 .0000 .000002 .0010 .0144 # .0040 
86 1000 50 .0000 .000002 .0025 .0237 .0049 
87 1000 50 .0000 .000002 .0100 .0522 .0050 
88 1000 50 .0000 .000005 .0010 .0144 .0074 
. 89 1000 50 .0000 .000005 .0025 .0237  .0109 
90 1900 50 .0000 .000005 .0100 .0522 .0121 
91 1000 50 .0050 .000001 .0010 .0144 .0066 
7 92 106. 50 .0050 .000001 .0025 .0237  .0074 
93 1000 50 .0050 .000001 .0100 .0522 .0075 
94 1000 50 .0050 .000002 .0010 .0144 .0079 
95 1000 50 .0050 .000002 .0025 .0237  .0095 
96 1000 50 .0050 .000002 .0100 .0522 .0099 
97 1000 50 .0050 .000005 .0010 .0144 .0105 
98 1000 50 .0050 .000005 .0025 .0237 .0149 
99 1000 50 .0050 .000005 .0100 .0522 .0168 
100 81000 50 .0075 .000001 .0010 .0144 .0088 
101 1000 50 .0075 .000001 .0025 .0237  .0098 
102 ~=1000 50 .0075 .000001 .0100 .0522 #.0100 
103. =1000 50 .0075 .000002 .0010 .0144 .0099 
104 1000 50 =.0075 .000002 .0025 .0237  .0118 
105 =1000 50 =©.0075 = .000002 .0100 .0522 .0124 
106 =©1000 50 .0075 .000005 .0010 .0144 .0122 
107. = 1000 50 .0075 .000005 .0025 .0237  .0168 
108 1000 50 .0075 .000005 .0100 .0522 .0190 
109 =1000 100 .0000 .000001 .0010 .0105 .0024 
110 i000 100 =..0000 §=.000001 .0025 .0175 .0024 
121 1000 100 =.0000 .000001 .0100 .0398 .0025 
. 112 ~=1000 100 =.0000 .000002 .0010 .0105 .0046 
113. 1000 100 =©.0000 = .000002 .0025 .0175 .0046 
114 =1000 100 =—.0000)3=—. 000002) «.0100=—Ss «. 0398 )~=—. 0050 
115 = 1000 100 .0000 .000005 .0010 .0105 .0086 


116 =1000 100 »=.0000 .000005 .0025 .0175 .0086 
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DEFINITION OF PARAMETERS 


SAMPLING FREQUENCY 


SAMPLE SIZE 


INHERENT PROBABILITY OF THE PROCESS PRODUCING A DEFECTIVE ITEM 


RATE AT WHICH THE PROBABILITY OF PRODUCING A DEFECTIVE 


ITEM INCREASES OVER TIME 
= DESIRED CENTERLINE VALUE ON CONTROL CHART (ROUGH ESTIMATE OF 


TiiE PROCESS AVERAGE FRACTION DEFECTIVE) 
COMPUTED UPPER CONTROL LIMIT FOR CONTROL CHART 


RESULTING ESTIMATE OF THE PROCESS AVERAGE FRACTION DEFECTIVE 


. 000005 


- 000001 
. 009001 


-& 9001 
. 0 0002 
. 000002 
. 000002 
. 000005 
. 000005 
. 000005 
. 900001 
. 000001 
. 000001 
. 000002 
. 000002 
. 000002 
. 000005 
. 000005 
. 000005 
. 000001 
. 000001 
. 000001 
. 000002 
- 000002 
. 000002 
. 000005 
. 000005 
. 000005 
. 000001 
. 000001 
. 000001 
. 000002 
. 000002 
. 000002 
. 000005 
. 000005 
- 000005 
- 000001 
. 000001 
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DEFINITION OF PARAMETERS 


SAMPLING FREQUENCY 

SAMPLE SIZE 

INHERENT PROBABILITY OF THE PROCESS PRODUCING A DEFECTIVE ITEM 
RATE AT WHICH THE PROBABII-ITY OF PRODUCING A DEFECTIVE 

ITEM INCREASES OVER TIME 

DESIRED CENTERLINE VALUE ON CONTROL CHART (ROUGH ESTIMATE OF 
THE PROCESS AVERAGE FRACTION DEFECTIVE) 

COMPUTED UPPER CONTROL LIMIT FOR CONTROL CHART 

RESULTING ESTIMATE OF THE PROCESS AVERAGE FRACTION DEFECTIVE 


m s a b cl UCL P 
1000 150 = .0075. —S— . 000001 =.0100 = .0344 = .0100 
1000 150 =©.0075 .000002 .0010 .0087 # .0096 
1000 150 .0075 .000002 .0025 .0147 .0110 
1000 150 = .0075 =. 000002.—S «0100S «© 0344 ~—(iw 012 
1000 180 «.0075—Ss «000005 =.0010 .0087 = .0113 
1000 150 =©=.0075 .000005 .0025 .0147 #.0137 
1000 150 =©=.0075 = .000005 .0100 .0344 .0189 
2000 50 =. 0000-—S«. 000001 =. 0010S. 0144 ~—«. 0040 
2000 50 .0000 .000001 .0025 .0237 + .0049 
2000 50 .0000 .000001 .0100 .0522 .0050 
2000 50 =©=.0000 .000002 .0010 .0144 # .0065 
2000 50 .0000 .000002 .0025 .0237 .0091 
2000 50 .0000 .000002 .0100 .0522 .0098 
2000 50 .0000 .000005 .0010 .0144 .0108 
2000 50 .0000 .000005 .0025 .0237 .0172 
2000 50 .0000 .000005 .0100 .0522 # .0217 
2000 50 .0050 .000001 .0010 .0144 #.0079 
2000 50 .0050 .000001 .0025 .0237 #.0095 
2000 5@ .0050 .000001 .0100 .0522 # .0099 
2000 59 . 0950 - 000002 . 0010 ~ 0144 . 0097 
2000 50 .0050 .000002 .0025 .0237 = .0133 
2000 50 .0050 .000002 .0100 .0522 # .0146 
2000 50 .0050 .000005 .0010 .0144 .0134 
2000 50 =6©.0050 =.000005 .0025 .9237 .0203 
2000 50 .0050 .000005 .0100 .0522 .0254 
2000 50 .0075 .000001 .0010 .0144 .0099 
2000 50 .0075 .000001 .0025 .0237 .0118 
2000 50 0075 - 000001 . 0100 . 0522 0124 
2000 50 .0075 .000002 .0010 .0144 .0115 
2000 50 .0075 .000002 .0025 .0237 0153 
2000 50 .0075 .000002 .0100 .0522 .0169 
2000 50 .0075 .000005 .0010 .0144 .0149 
2000 50 .0075 .000005 .0025 .0237 + .0218 
2000 50 .0075 .000005 .0100 .0522 .0272 
2000 100 =.0000 .000001 .0010 .0105 .0046 
2000 100 .0000 .000001 .0025 .0175 .0046 
2000 100 .0000 .000001 .0100 .0398 .0050 
2000 100 .0000 .000002 .0010 .0105 = .0076 


2000 100 =.0000 .000002 .0025 .0175 = .0076 











DEFINITION OF PARAMETERS 


m = SAMPLING FREQUENCY 

s = SAMPLE SIZE 

a = INHERENT PROBABILITY OF THE PROCESS PRODUCING A DEFECTIVE ITEM 

b = RATE AT WHICH THE PROBABILITY OF PRODUCING A DEFECTIVE 
ITEM INCREASES OVER TIME 

cl = DESIRED CENTERLINE VALUE ON CONTROL CHART (ROUGH ESTIMATE OF 
THE PROCESS AVERAGE FRACTION DEFECTIVE) 

UCL = COMPUTED UPPER CONTROL LIMIT FOR CONTROL CHART 


p = RESULTING ESTIMATE OF THE PROCESS AVERAGE FRACTION DEFECTIVE 


-000002 .0100 .0398 .0097 
.000005 .0010 .0105 .0117 
.000005 .0025 .0175 .0117 
-000005 .0100 .0398 .0192 
.000001 .0010 .0105 .0086 
.000001 .0025 .0175 .0086 
.000001 .0100 .0398 .0099 
.000002 .0010 .0195 .0109 
.000002 .0025 .0175 .0)09 
.600002 .0100 .0398 .0143 
.000005 .0010 .0105 .0145 
.000005 .0025 .0175 .0145 
.000005 .0100 .0398 .0222 
.000001 .0010 .0105 .0106 
.000001 .0025 .0175 .0106 
.000001 .0100 .0398 .0123 
-000002 .0010 .0105 .0125 
.000002 .0025 .0175 .0125 
.000002 .0100 .0398 .0165 
.000005 .0010 .0105 .0159 
. 000005 .0025 .0175 .0159 
- 000005 .0100 =.0398 8.0237 
. 000001 9010 .0087 .0042 
.000001 .0025 .G147 .0048 
.000001 .0100 .0344 .0050 
.000002 .0010 .0087 .0063 
.000002 .0025 .0147 .0082 
.000002 .0100 .0344 .0099 
.000005 .0010 .0087 .0093 
.000005 .0025 .0147  .0121 
.000005 .0100 .0344 .0201 
.000001 .0010 .0087 .0078 
.000001 .0025 .0147 = .0090 
.000001 .0100 .0344 .0100 
.000002 .0010 .0087 .0092 
.900002 .0025 .0147 .0114 
.000002 .0100 .0344 ~~ .0147 
.900005 .0010 .0087 .0120 
.0000C5 .0025 .0147 .0149 
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DEFINITION OF PARAMETERS 





m = SAMPLING FREQUENCY 
s = SAMPLE SIZE 
a = INHERENT PROBABILITY OF THE PROCESS PRODUCING A DEFECTIVE ITEM 
b = RATE AT WHICH THE PROBABILITY OF PRODUCING A DEFECTIVE 
ITEM iNCREASES OVER TIME 
cl = DESIRED CENTERLINE VALUE ON CONTROL CHART (ROUGH ESTIMATE OF 
THE PROCESS AVERAGE FRACTION DEFECTIVE) 
UCL = COMPUTED UPPER CONTROL LIMIT FOR CONTROL CHART 
p = RESULTING ESTIMATE OF THE PROCESS AVERAGE FRACTION DEFECTIVE 
Case m Ss a b cl UCL p 
234 2000 150 =©.0050 .000005 .0100 .0344 .0231 
235 2000 150 = 60075. =. 000001 =.0010 =.0087 =. 0096 
236 = §=2000 150 = .0075 .000001 .0025 .0147 + .0110 
237 ~=2000 150 =©.0075 .000001 .0100 .0344 .0124 
238 2000 150 .0075 .000002 .0010 .0087 .0108 
239 2000 150 .0075 .000002 .0025 .0147 #.0130 
240 =2000 150 =©=.0075 =.000002 .0100 .0344 #.0169 
241 2000 150 .0075 .000005 .0010 .0087 + .0133 
242 2000 150 =©.0075 .000005 .0025 .0147 # .0162 
243 =. 2000 150 =©.0075 .000005 .0100 .0344 .0246 
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